Negative group velocity and three-wave mixing in 

dielectric crystals 



Vitaly V. Slabko, 1 Sergey A. Myslivets, 2 Mikhail I. Shalaev, 1 and Alexander K. 

Popov 3 '* 

1 Siberian Federal University, 660041 Krasnoyarsk, Russia 
2 Institute of Physics of the Russian Academy of Sciences, 660036 Krasnoyarsk, Russia 
3 University of Wisconsin-Stevens Point, Stevens Point, WI 54481, USA 
* Corresponding author: apopov@uwsp.edu 



Extraordinary features of optical parametric amplification of Stokes electro- 
magnetic waves are investigated, which originate from the three-wave mixing 
of two ordinary electromagnetic and one backward phonon wave with negative 
group velocity. A similarity with the counterpart in the negative-index plas- 
monic metamaterials and differences with those utilizing contra-propagating 
ordinary electromagnetic waves as well as electromagnetic and acoustic 
phonon waves are shown. They stem from the backwardness of the optical 
phonons in crystals with negative dispersion. Nonlinear-optical photonic 
devices with the properties similar to those predicted for negative-index 
metamaterials are proposed. ©2013 Optical Society of America 

OCIS codes: 190.4223, 290.5910, 260.1180, 350.3618, 230.4320. 
1. Introduction 

Optical negative-index (NI) metamaterials (NIMs) form a novel class of electromagnetic 
media that promises revolutionary breakthroughs in photonics pQ. The possibilities of such 
breakthroughs originate from backwardness, the extraordinary property that electromagnetic 
waves acquire in NIMs. Unlike ordinary positive-index materials, the energy flow, S, and 
the wave-vector, k, become counter-directed in NIMs, which determines their extraordinary 
linear and nonlinear optical (NLO) propagation properties. Usually, NIMs are nanostructured 
metal-insulator composites with a special design of their building blocks at the nanoscale 
[2jl3] . Generally, the metal component imposes strong absorption of optical radiation in NIMs, 
which presents a major obstacle towards their numerous exciting prospective applications. 
Extraordinary features of coherent NLO energy conversion processes in NIMs that stem 
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from wave-mixing of ordinary and backward electromagnetic waves and the possibilities 
to apply them for compensating the outlined losses have been shown in [HE]. Essentially 
different properties of three- wave mixing (TWM) and second harmonic generation have been 
revealed in [6l[7] . Herein, we propose and investigate a different scheme of TWM of ordinary 
and backward waves (BW). It builds on the stimulated Raman scattering (SRS) where two 
ordinary electromagnetic waves excite backward elastic vibrational wave in a crystal, which 
results in TWM. The possibility of such BWs was predicted by L. I. Mandelstam in 1945 [5], 
who also had pointed out that negative refraction is a general property of the BWs. This letter 
is to show the possibility to replace negative index composites with the extensively studied 
ordinary crystals and thus to simulate unparallel properties of coherent NLO energy exchange 
between the ordinary and backward waves. We show extraordinary nonlinear propagation 
and output properties of the Stokes electromagnetic wave in one of two different coupling 
geometries, both utilizing backward elastic waves. The possibility to localize coherent energy 
conversion and to fit it in to the crystal of given thickness is shown. Such unusual properties 
are in a striking contrast with those attributed to the counterparts in the standard schemes 
that build on the coupling of co-propagating photons and phonons j9j[]Il|. They are also 
different from the properties of the phase-matched mixing of optical and acoustic waves for 
the case when the latter has energy flux and wave vector directed against those of one of the 
optical waves [TT]. The revealed properties can be utilized for creation of optical switches, 
filters, amplifiers and cavity-free optical parametric oscillators based on ordinary Raman 
crystals without the requirement of periodically poling at the nanoscale [12] (and references 
therein) . 

2. Optical phonons, negative group velocity and phase matching options 

The dispersion curve u(k) of phonons in the crystals containing more than one atom per 
unit cell has two branches: acoustic and optical. For the optical branch, the dispersion is 
negative, in the range from zero to the boundary of the first Brillouin's zone (Fig. [T]) and the 
group velocity of optical phonons, v-f", is antiparallel with respect to its wave- vector, k£ , 
and phase velocity, v^ 1 , because v gr = du(k)/dk < 0. Optical vibrations can be excited by 
the light waves due to the two-photon (Raman) scattering. The latter gives the ground to 
consider crystal as the analog of the medium with negative refractive index at the phonon 
frequency and to examine the processes of parametric interaction of the three waves, two 
of which are ordinary electromagnetic waves and the third is the wave of elastic vibrations 
with the directions of the energy flow and of the wave-vector opposite to each other. 
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Fig. 1. Negative dispersion of optical phonons and two phase matching options: (a) - co- 
propagating, (b) - contra-propagating fundamental, control, and Stokes, signal, waves. Insets: 
relative directions of the energy flows and the wave- vectors. 

3. Basic equations 

We will examine only lowest-order Raman process [9j[T0]. The interacting waves are given 
by the equations 

E lt , = (l/2)e l ,.(z,t)e ikt " s - iUl " t + c.c., 
Q v = (l/2)Q(z,t)e lkvZ - iuJvt + c.c. 

Here, Ej )g , Q, loi >s>v and ki >8iV are the amplitudes, frequencies and wave-vectors of the funda- 
mental, Stokes and vibrational waves; Q v (z, t) = y/px(z, t); x is displacement of the vibrating 
particles, p is the medium density. With account for the energy and momentum conservation, 

ui = u s + u v (k v ) , ki = k s (u s ) + k v , 
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one obtains the following equations for the slowly varying amplitudes in the approximation 
of the of first order of Q in the polarization expansion: 

dz + vf dt l kc 2 ■ dQ £a ^ 
de s 1 de s nu 2 s AJ da 

^ + + X = i—N^e l£ *. (2) 

dz Vv dt Tvf, r 4:U v Vv T dQ s 

Here, vf r s v are the projections of the group velocities of the fundamental, Stokes and vibration 
waves on the z-axis, N is the number density of the vibrating molecules, a is the molecule 
polarizability, r is phonon lifetime, uq is phonon frequency for k v = 0. The dispersion u v (k v ) 
can be approximated as [9] 

Then, in the vicinity of k v = 0, velocity v 9 J is given by: 

< r = -(3k v /u v = -[3/v p v h , 

where v% h is the projection of the phase velocity of the vibrational wave on the z-axis and /3 
is the dispersion parameter for the given crystal. 

For the sake of clarity, the continuous wave case and the approximation of the constant 
field Ei is considered. The latter is appropriate for the relatively week Stokes and vibrational 
waves. Then Eqs. ([I])-© take the form: 

dQ/dz = -ig\£* s - Q/(Tv 9 J),de s /dz = ig 2 Q*, 

where g% = —N(da/dQ)ei/(4u v v 9 , r ) and g 2 = ( r nuj 2 s /k s c 2 )N{da/dQ)ei. In the case of 
Fig. Ufa), Eqs. ()3]) exhibit three fundamental differences as compared with TWM of co- 
propagating waves in ordinary materials: an opposite sign with gi which stems from v 9 J < 0, 
an opposite sign with Q/(tv 9T ) because the phonon flow is against the z-axis, and the bound- 
ary conditions for Q to be defined at z — L, i.e. at the opposite edge of the slab as compared 
to that for e s . This leads to fundamental changes in their solutions and, consequently, in 
the spatial and output behavior of the Stokes signal. Alternatively, in the given constant 
ei approximation, the equations become identical and the behavior standard for the case of 
Fig. Lib). 

4. Extraordinary properties of coherent energy conversion from fundamental 
to Stokes electromagnetic waves 

The solution to Eqs. ([3]) is found in the form: 

Q* = A ie ^ z ' + A 2 e^ z ', e s = A 3 e^ z ' + A A e^ z ' , (3) 
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where /?i j2 = 1 =F iR, R = \/ 9i92^ — 1, -2' = ^Ap, £p = — 2v% r T. The amplitudes Ai_4 
and their relationships are determined by the boundary conditions. Transmission factors for 
co-propagating, Tj^(z), and counter-propagating (g 2 < 0), TJ^(z), fundamental and Stokes 
waves are found as 

2 

(4) 



e z ' {R cos [R {V - z')} + sin [R (V - z')} } 



Rcos(RL') + sm (RU) 



where L' = L/l p , = \e s (z)/e s (z = 0)| 2 and = \e s (z)/e s (z = L)\ 2 . Equations (j4j) and 
§5§ display spatial distributions which are controlled by the field S[ and are in a strict contrast 
to each other. It is explicitly seen for the ultimate loss- free case (l p —> 00). Then 

Tf{z = L) l/cos 2 (^L), (6) 
Tj\z = 0) [exp(2^|L)]/4, (7) 

where g = y/g~lg~2. Equation (Q depicts a series of sharp giant resonance enhancements of 
the output signal for g (2j + l)n/2L, (j = 0, 1, 2...). On the contrary, the coupling scheme 
of Fig. QJb) is equivalent to scattering on acoustic phonons and on optical phonons with 
positive group velocity. Correspondingly, Eq. (J7|) displays typical exponential growth with 
no resonances with respect to intensity of the fundamental control field. In general case, the 
denominator in Eq. (j3J) can be turned to zero if g 2 ti > 1. The threshold value of intensity of 
the control field is 

/rnin = {c7l s \ s0 tU v / 8iT%t) \N da / dQ\~ 2 , (8) 

where \ s q is Stokes wavelength in the vacuum. For a given intensity of the control field 
h > Imini the crystal thickness corresponding to the first resonance is 

L'=[n- tan" 1 (R)] / R. 

Figure EJ^a) depicts transmission in the vicinity of the first resonance and Fig. EJ^b) - its 
position. In the resonance, T}} — > oo, which is due to the approximation of constant control 
field. Conversion of the control field to the Stokes one and excited molecule vibrations would 
lead to saturation of the control field and thus limits the maximum achievable amplification. 
Strong amplification in the maximums indicates the possibility of self-oscillations and thus 
creation of mirrorless optical parametrical oscillator with unparalleled properties. 

In the case of Fig. [TJb), k s < and the denominator in Eq. cannot equal to zero. 
As outlined, this results in exponential spatial dependence with no resonances depicted in 
Fig. E^c). Figure E(d) shows that, in the vicinity of the resonance, three-wave coupling of 
waves with co-directed wave vectors and contra-directed energy flows of vibrational and 
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Fig. 2. (a) Transmission of the Stokes wave Tf(z = L) vs intensity of the fundamental control 
field in the vicinity of the first resonance (co-propagating geometry), (b) The crystal thickness 
corresponding to the first output resonance for the Stokes wave vs intensity of the control field 
for the co-propagating geometry; J min is the threshold intensity, (c) Transmission Tj^(z = 0) 
vs intensity of the fundamental beam (contra-propagating geometry), (d) Comparison of the 
output intensities of the Stokes wave vs intensity of the control field for co- and contra- 
propagating coupling geometries. 



6 



Stokes waves provides for essentially higher efficiency of coherent energy conversion than in 
the standard schemes. Figures W[b, d) indicate the possibility to fit in the effective conversion 
length within the crystal of a given thickness and to significantly concentrate the generated 
Stokes field nearby its output facet. 

Estimations made for the model, which is characteristic for the diamond crystal, u v = 1332 
cm -1 and vibrational transition width (cr) _1 = 1.56 cm -1 [TBTnS] . show that the required 
excitation intensities are above the typical crystal breakdown threshold in the continuous 
wave regime. However, the cross-section of the Raman scattering is inverse proportional 
to the squared frequency offset from the intermediate single-photon resonance. Therefore, 
the threshold intensity l m in given by Eq. (JHJ) can be reduced by seven to eight orders com- 
pared to its off- resonance values by approaching such a resonance. The required intensity is 
also achievable with the commercially available pulsed lasers, especially for the pulse length 
shorter than the phonon relaxation rate. 

5. Conclusions 

To conclude, we propose extraordinary coherent energy conversion process attributed to 
three-wave mixing of ordinary and backward waves. It rests on the excitation of optical 
phonons with negative dispersion and can be implemented in the routinely fabricated crystals 
extensively studied by the means of Raman spectroscopy. The investigated features can be 
employed for creation of a family of the unique nonlinear-optical photonic devices, such as 
frequency-tunable switches, filters, amplifiers and miniature mirrorless optical parametric 
oscillators. Earlier, such a possibility was proposed based on negative index metamaterials 
which are much more difficult to fabricate and are inherently strongly lossy. The revealed 
features stem fundamentally from the opposite directions of the phase velocity and the energy 
flow, which is the characteristic property of optical phonons having negative dispersion, 
dco(k)/dk < 0. 
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